Three-Body Halo Fragmentation: Polarization effects 



E. Garrido, D.V. Fedorov and A.S. Jensen 

Institute of Physics and Astronomy, 
Aarhus University, DK-8000 Aarhus C, Denmark 
(February 9, 2008) 

Momentum distributions of particles from nuclear break-up of fast three-body halos are calcu- 
lated for polarized beam and outgoing fragments. These momentum distributions are combined 
in observable quantities which emphasize the two-body correlations in the initial three-body struc- 
ture as well as in the final two-body system. Applying the formalism to 11 Li we find a sensitive 
, dependence on the p-wave content and on the position of the p-resonances in 10 Li. Polarization 

^\ • experiments are then a new tool to study correlations in nuclear halos. 

G\ ; 

PACS. 25.60.+V - Reactions induced by unstable nuclei. 
bJ[)| PACS. 24.70.-fs - Polarization phenomena in reactions. 

£3 ! PACS. 21.45. +v - Few-body systems. 
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Introduction. The unusual halo nuclei appear under conditions realized at the neutron dripline The structure 
of these nuclei is dominated by a configuration where an ordinary nucleus (core) is surrounded by weakly bound 
neutrons. For Borromean systems J^||, three particles without bound two-body subsystems, specific correlations 
must be present to provide the total binding. The nature of these correlations is rather difficult to investigate and 

■ the data have almost always been consistent with uncorrelated motion of particles around the core. The mixture of 
£f~) \ s- and p-states in the neutron-core subsystem is an obvious manifestation of the correlations. However its existence 

■ has been rather difficult to confirm experimentally. 

The major source of detailed information is the measurements of particle momentum distributions after fragmen- 
tation reactions of these halo nuclei . In the typical experimental setup two of the three constituent particles are 
detected after the high energy halo fragmentation in anticoincidence with the third particle, which has been instantly 
captured by the target. The sudden approximation or the Serber mechanism, equivalent to instantaneous cease of the 
interactions between the constituent particles, has been proved to be adequate in data interpretation Further- 
1 , ■ more, it has been pointed out recently that inclusion of the interaction between the two remaining particles greatly 
improves the model 

A polarized beam interacting with a target is a tool sensitive to components of the wave function with non-zero 
relative angular momentum. In particular for 11 Li, the resulting particle momentum distribution should be sensitive 
to the amount of p-state admixture in the neutron-core relative state. More details and strong signals are expected 
when the polarization of the fragments is simultaneously measured. 

The purpose of this letter is to investigate the effects of polarization on halo nuclear reactions in view of extracting 
information about the correlations in these nuclei. We shall first briefly sketch the new method for computation 
of fragment momentum distributions arising from polarized beam experiments. We shall investigate the effects of 
polarized beams and unpolarized fragments, as well as polarized beams and polarized fragments. 

Method. In the sudden approximation the transition matrix is proportional to the overlap between the initial and 
the final wave functions. This approximation is suitable for neutron halo fragmentation, as first suggested for the 
deuteron stripping by Serber fll|| and later used in many connections [fl2f . The projectile wave function is denoted 
by $ JM (x, y), where J and M are the total spin and its projection. The coordinates x and y are the usual Jacobi 
coordinates |2[|, where x is drawn between the two particles surviving after the fragmentation. The final state wave 
function is given by e l ^ y ' y Xsl J2 Sx a x ( s li s 2; &i, o-2\s x a x )w Sx<7x (k x , x), where w s *' J * (k x , x) describes the final two-body 
system with spin s x and projection o~ x . The quantum numbers s\, S2, y , ci, <J 2 , and a y are the spins and projections 
of the three particles, and the plane wave describes the motion of the center of mass of the two-body system. The 
plane waves || correctly describe the shapes of the high energy momentum distributions in the region of the first 
maximum. By squaring this overlap, and subsequently summing up over the initial and final magnetic substates we 
get the following expression for the cross section or momentum distribution 
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where tVjnit (M) , Wa n (ai), and W& a (<J2) give the probability of occupation of each magnetic substate in the initial 
and final states. When the beam is not polarized and there is no detection of final spin states these probabilities take 
the values 1/(2J+1), l/(2si+l), and l/(2s2+l), and we obtain the unpolarized expression. 

From x and y we define the hyperspherical coordinates p = y 'x 2 + y 2 , a = arctan(x/y), and fl x and O y giving the 



directions of x and y. Analogously, in momentum space we define k = yk 2 + k 2 , a K = axctan(k x / k y ) , fifcj,, an d ^fc„- 
The three-body wave function is obtained by solving the Faddeev equations in coordinate space Ei] . This function 
is then expanded in terms of the complete set of hyperspherical harmonics Y^"'^ (a, f2 x , Fl y ), where l x and i y are the 
orbital angular momenta associated to the Jacobi coordinates x and y, and L is the coupling of both momenta. The 
quantum number K is usually called the hypermoment. The function $! JM (x, y) can then be expanded as 



* JM (x, y) = ]T F«JX (P) («' Vy) ® xLs y (2) 
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where Xs x Sy is the spin function obtained by coupling the two-body spin s x and the spin s y of the third particle to 
the total spin S. 

We expand u> Sx<Tx (k^, x) in partial waves 

-7— E u3 £s x ( k x> x )yZe x sA n x) ^^xS x ]m tx a x \j x m x )i l x Y^ mix {VL k:c ), (3) 
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where the radial functions u J f s (k x , x) are obtained by solving the Schrodingcr equation with the corresponding two- 
body potential. The quantum numbers contained in eq.(^) correspond to the ones in the two-body subsystem in 
eq.(||). Following this procedure we incorporate Final State Interactions (FSI) into the calculation. When FSI is not 
considered the function ■u; Sx<Tx (k^, x) reduces to a plane wave, and the overlaps in eq.(Q) become the Fourier transform 
of the three-body wave function. 

First we consider the case of 100% polarized beam {W- m {t{M') = $m'm)- After substituting eqs.(||) and (g) into 
eq.(Q), and integrating analytically over VL ky and ipk x {^k x = {@k x , yfc x ))j we get the expression 
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where k y — \k y \, k x — \k x \ sin 6*fc x , fe| = |k x | cos 9k x , and a = y/2a + 1. The functions (k, a K ) are only available 

as numerical functions. By integration over k y and either fc 1 - or kx, we obtain respectively Qf(fcl) or Q.% 1 (k^r), the 
one-dimensional or two-dimensional momentum distribution (note that in our description there is no preferred axis, 
and therefore the one-dimensional longitudinal and transverse momentum distributions are identical). 

From the definitions of the coordinates it is clear that Q* f 2 expresses the relative momentum distribution of the 
two particles in the final state. To compare with the experimental data one has to refer the momentum distributions 
to the center of mass of the three-body system. The appropriate expressions, used in the numerical calculations, are 
rather complicated but the qualitative features remain unchanged. Details about these expressions will be given in a 
forthcoming paper. 

The summation over M in eq.(^) leads to the unpolarized momentum distribution, that amounts to keeping only 
the I — term. The terms with 1^0 appear only when the beam is polarized, and only £ x =^ waves contribute. 
We now define the asymmetry 



A f M ' = ^ nL ' * = 1 ' 2 ( 5 ) 

that vanishes when the beam is not polarized. The deviations from zero are then produced by the I ^ terms in 
Qf 1 , or, in other words, by the t x ^ waves contained in the neutron-core subsystem. 
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Finally, together with a polarized beam, we consider the possibility of measuring the polarization of one of the 
fragments in the final state. Then the momentum distribution is similar to eq.(|j), but containing additional geometrical 
factors involving the z-projection a of the spin of the polarized particle in the final state (the expressions are rather 
extended, and the details will be published in a coming paper). From the one and two-dimensional momentum 
distributions Qf Icr (kx) and Q^^^h^r) we then introduce the asymmetry 
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where again the deviations from zero are produced by the polarizations. However, in this case also s-waves contribute. 

Results. We apply the method to the three-body halo 11 Li ( 9 Li+n+n) fragmentation, where one of the neutrons 
is suddenly captured by the target. The two-body potentials used in the calculations are given in ref. p3[ , where 
the nucleon-nucleon potential is fitted to low energy s- and p-wave nucleon-nucleon scattering data, and the neutron- 
core potential is adjusted to reproduce the experimental binding energy and mean square radius of the three-body 
system. In all the calculations the lowest virtual s state in 10 Li is at 50 keV, and the lowest p-resonance is at 0.5 
MeV, according to the experimental data ||. The rest of the resonances occurs at different energies, modifying the 
p-content in the two-body subsystem and keeping the right three-body binding energy. 
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FIG. 1. Upper part: Asymmetry A 2 2 for neutron momentum distributions from lx Li fragmentation reaction with polarized 
beam. Three different contents of p-wave in the neutron-core system are considered: 35% (solid), 18% (long dashed), and 4% 
(dashed). The curve for 35% of p-wave without final state interactions (dot-dashed) is also shown. Lower part: The same as 
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in the upper part for the asymmetry 2 ' 2 2 . Inner part: Unpolarized one-dimensional neutron momentum distribution from 
neutron removal process in 1 Li fragmentation. The neutron momentum p is relative to the center of mass of the three-body 
system. 
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Let us first focus on the asymmetry A\ 2 for momentum distributions of the surviving neutron. This asymmetry 
is shown in the upper part of fig. I, where three curves corresponding to different p-wave content in the neutron- 9 Li 
channel (4, 18, and 35%) are shown. 

From this part of the figure we can conclude the following: 

(i) The deviations from zero of the asymmetry are produced by the p-waves in the neutron- 9 Li subsystem. A 
significant dependence on the p-wave content of the two-body wave function is observed. Since the asymmetry is 
zero when no p-waves are present, one could expect smaller oscillations for small p-wave content. However, p-waves 
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are also contained in the unpolarized momentum distribution, in such a way that the less the p-wave content, the 
narrower the momentum distribution (see the inner part of fig.l). Therefore the denominator of the asymmetry tends 
to increase the amplitude of the oscillations for low p-wave content in the two-body wave function. It is then possible, 
as in fig.l, to have oscillations with similar amplitude for different values of the p-wave content. 

(ii) If the two-body system has a resonance at a given energy E res , then the enhanced overlaps in eq.(|l|) produce 
a bump in the momentum distribution at k = \/2fiE ICS (p is the reduced mass of the two particles). For the same 
reason the difference between Qf 1 and the unpolarized momentum distribution is expected to be larger in the regions 
close to the resonances in the two-body system. Therefore the asymmetry should present larger oscillations for the 
momentum values corresponding to p-resonances in 10 Li. In particular, the p-resonance at 0.5 MeV corresponds to a 
momentum of 30 MeV/c, that matches with the position of the first peak in the asymmetries. The second p-resonance 
changes its position from case to case, being placed at 1.5 MeV (k ~ 50 MeV/c) for the solid curve, 2 MeV (k ~ 60 
MeV/c) for the long-dashed curve, and 6 MeV (k ~ 100 MeV/c) for the dashed line. From the figure it is then clear 
that the asymmetry gives an indication about the level structure of 10 Li. 

(iii) If we consider the three particles in the initial state coupled to total orbital angular momentum L = 0, the 
polarization of the beam affects only the spin part of the wave function. If then none of the final spin projections are 
measured, neglecting the interaction between 9 Li and the remaining neutron (no FSI) leads to momentum distributions 
that are equal to the unpolarized ones, and the asymmetry therefore vanishes (this can be seen analytically by 
inserting eqs.(||) and into eq.(]l|) for L = 0). In our case, the n Li wave function contains more than 95% of L = 
configuration, and therefore the asymmetry neglecting FSI is almost zero (see the dot-dashed line in the figure). We 
can therefore conclude that FSI is necessary for the correct description of the asymmetry. 




We turn now to the case where the polarization of the surviving core also is measured. We assume that the core 
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is longitudinally polarized and consider the asymmetry A\ 2 ' 2 2 . In the lower part of fig.l we plot this asymmetry 
for the same three cases we considered above. The three points we already emphasized are valid again. First a clear 
dependence on the p-wave content is observed even though the s-waves also appear in the numerator of the asymmetry. 
The effect of s-waves is however confined to the low momentum region, and for momenta larger than 20 MeV/c only 
p-waves contribute. Therefore the minima and maxima of the asymmetry again match with the p-resonances in 10 Li. 
Finally the dot-dashed line shows that FSI is also essential in this case. However as the polarization of the core is 
now measured, even L — terms in the three-body wave function contribute to the asymmetry. This asymmetry is 
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then an observable with similar characteristics as A\ 2 , but with the advantage that the oscillations are larger, and 
therefore easier to observe experimentally. 

In fig. 2 we show the same calculations as in fig.l, but for two-dimensional neutron momentum distributions. The 
curves show features similar to the curves in fig.l. Our conclusions about the behaviour of the asymmetries are 
therefore rather general and not restricted to only one kind of momentum distributions. 

Summary and conclusions. The aim of this letter is to investigate polarizations in three-body halo fragmentation 
reactions at high energies and light targets. We use the sudden approximation where the initial state is a genuine 
three-body wave function and the final state is a distorted wave in the same potential between the two surviving 
particles as used in the three-body calculation. We have defined several asymmetries from which the polarization 
effects can be analyzed. 

We applied the method to 11 Li fragmentation, and computed one and two-dimensional neutron momentum distri- 
butions. We have found that the polarization observables are more sensitive to the p-wave content in the neutron-core 
subsystem than unpolarized cross sections. Furthermore they provide a unique information of the level structure of 
the two-body subsystem. We also observed that FSI is the main cause for the behaviour of the asymmetries, and 
FSI must necessarily be included in the calculation. These features have been found for fragmentation reactions with 
polarized beam, and also when simultaneously the polarization of the core is measured. The advantage of the second 
case is that the amplitude of the oscillations observed in the asymmetries is clearly larger. Finally we have shown 
that these conclusions are valid for both one and two-dimensional momentum distributions. 
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